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Abstract 

We perform constant pressure Monte Carlo simulations of a spin-analogous model which describes cou- 
pled spatial and magnetic degrees of freedom on an fee lattice. Our calculations qualitatively reproduce 
magnetovolume effects observed in some rare earth manganese compounds, especially in the anti-Invar 
material YMn2. These are a sudden collapse of the magnetic moment which is connected with a huge 
volume change, and a largely enhanced thermal expansion coefficient. 

INTRODUCTION 

The observation that some transition metal alloys and complexes have an anomalously low or high 
thermal expansion coefficient is widely referred to as the Invar or anti-Invar effect, respectively. Since 
the discovery of Fe-Ni Invar one century ago, uncounted investigations in this topic have taken place 
(for recent reviews see [Q). A major advance has been achieved by ab initio band structure calculations 
(e. g. d, ^). KKR-CPA calculations of FegsNias reveal two nearly degenerate minima on the binding 
surface, i.e. the energy as a function of lattice constant and mean magnetic moment. One at a large 
lattice constant and a large magnetic moment (high moment or HM state) and another at a smaller 
lattice constant and a smaller or vanishing magnetic moment (low moment or LM state). In this case, 
thermal excitation of the LM state can lead to a compensation of the volume expansion of the material. 
This scenario reminds of the phenomenological 2-7-states model Q introduced by Weiss in order to 
explain Fe-Ni Invar. However, most ab initio calculations are performed for zero temperature only. So 
the comparison with experimental data is hampered by the fact that the extrapolation of these results 
to finite temperatures is connected with severe methodological problems. So far mainly continuous 
Ginzburg Landau like spin models in combination with a Gaussian fluctuation theory have been used to 
achieve this goal j|] . Another approach is to use a Weiss-type model with its parameters adapted to ab 
initio results. This has been done for the classical Invar system Fe-Ni, and although we neglected the 
itinerant character of magnetism by using localized spins, we were able to reproduce qualitatively the 
major magnetovolume effects of this alloy 

However, it has still to be proven, whether this approach can be applied to other materials, especially 
to anti-Invar systems. One of the most prominent anti-Invar materials is the Laves phase compound 
YMn2. It occurs in the C15 structure with a frustrated antiferromagnetic (AF) spin order in the Mn 
sublattice. Around Tn ~ 100 K YMn2 undergoes a first order phase transition into a paramagnetic 
state, accompanied by a huge volume contraction of 5 % and a large thermal hysteresis of 30 K. Above 
Tn a large thermal expansion coefficient of 50 x 10~^ K^^ is encountered [Q. LMTO calculations reveal 
a nonmagnetic solution at a lower volume with a slightly higher energy than the AF ground state . 

Within the scope of this work, we will show that for a suitable choice of parameters the large abrupt 
volume change at Tn as well as the enhanced thermal expansion of YMn2 can be qualitatively reproduced 
by a Weiss type model assuming a HM-LM instability of the Mn atom. 
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THE MODEL 



Our model Hamiltonian consists of a magnetic part Hy^ describing the the magnetic properties of the 
system and a vibrational part responsible for the magnetovolume coupling: 

H = DY,Sf + jY^ SA. + U{nk,S,,Sk) . (1) 

i <i.k> <i,k> 



ifm is a spin-1 Ising Hamiltonian, also known as Blume-Capel model. We identify the spin states Si — 
with the atomic LM state, whereas Si = ±1 refers to the HM state. D denotes the crystal field, which 
imposes an energetic separation between HM and LM states; J is the magnetic exchange constant. The 
brackets indicate a summation over nearest neighbors. We did not consider a distance dependence of J so 
far, since we believe its influence to be small compared to the effect observed. Furthermore, compressible 
frustrated Ising antiferromagnets show additional effects like tetragonal lattice distortions §, that are 
not covered by experimental data in our case. The second part introduces pair interactions between 
the atoms, which depend on their distance r^^ and on their spin states. If both of the interacting atoms 
are in the HM state, we employ HM potentials that are characterized by a larger equilibrium distance 
than the LM potentials we use otherwise. So a system with a considerable concentration of LM atoms 
is supposed to have a smaller lattice constant than a pure HM system. For the sake of simplicity we 
choose Lennard- Jones interactions: 



SiSk — 

where eL.H denote the energy at the equilibrium neig hbor distances r^^ = h for LM and HM 

potentials, respectively. 



THE SIMULATION 



Our Monte Carlo routine consists of two local and one global update steps. For each lattice site we 
choose a new spin state Si G {0, ±1} using the Metropolis criterion. Afterwards a new trial position is 
elected out of a cube around the old position. Again, the new position is accepted with the probability 
max(l, exp(— /3Ai?)), where f3 is the inverse temperature. The size of the cube is given by the condition 
that about half of the propositions are to be accepted in order to improve convergence. When all 
atoms have been updated, the volume of the complete lattice is adapted by another Metropolis step, 
except that now the quantity AH = AH — N k^T In {V' /V) has to be considered, which accounts for 
the difference in translational entropy caused by the change of the volume from V to V (for more 
information on constant pressure Monte Carlo methods see |p^ ). Our calculations were performed on 
a 8'^ X 4 = 2048 site fee lattice with periodic boundary conditions at constant pressure P = 0. 80 000 
lattice sweeps were performed for each temperature. The first 40 000 lattice sweeps were used to allow 
the system to reach equilibrium, afterwards every tenth sweep data were gathered for the statistics. 

Although YMn2 occurs in the C15 structure, we chose instead an fee lattice in our simulations. This 
allows us to omit the Y species, which we do not expect to play an important role in the thermodynamics 
of the HM-LM transition. Furthermore we can easily relate our results to previous investigations Q. 
Since the fee Ising antiferromagnet has also a frustrated ground state spin structure as it is assumed 
for YMn2, simulating an fee structure is not a serious restriction. We used dn = 2.432 A and = 
2.397 A for the interatomic spacing between the Mn atoms in AF and paramagnetic YMn2. We chose 
El = ch = 30.86 mRy in order to achieve a reasonable elastic behavior. The exchange constant was 
set to ,/ = 0.459 mRy corresponding to a Neel temperature Tn = 126.5 K jl^. The crystal field D 
was varied in the range D = ... 3 J. We examined the energy per lattice site (e), the mean magnetic 
moment (S^) = {l/N "^iSf), the atomic volume {V) and the linear thermal expansion coefficient 
a = 1/(3V^) dV/dT (brackets denote thermal averaging). Additionally, for several values of D, we 
approximated ground state energy and average magnetic moment as a function of the lattice constant 
by exponentially cooling down the system to T = 4 K at constant volume. 
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RESULTS 



For D < 2 J, the ground state of the system is an AF HM state as can be seen from simple energetic 
considerations concerning Hamiltonian (^. Correspondingly, for D > 2 J, the ground state turns into 
a nonmagnetic LM state with a smaller volume. This is also verified by the plots of low temperature 
energy and magnetic moment vs. lattice constant shown in Fig. |^. 

Disallowing the LM spin state {D = — oo), around Tn — 127 K a first order phase transition occurs 
from a type 1 AF state to a paramagnetic state. This transition can be observed as a finite jump in 
the energy per lattice site (Fig. |^). Enlarging D allows thermal excitations of the 5* = spin state, 
leading initially to a continuous decrease of HM atoms with increasing temperature, which appears in 
our results for D — (Fig. ^). Due to the magnetovolume coupling introduced by H^, this is connected 
with a continuous decrease of the volume and, correspondingly, a reduction of the thermal expansion 
coefficient, as depicted by Fig. ^ The AF phase transition is hardly affected by the dilution of the 
magnetic system and remains of first order; only Tn is slightly reduced. We now find an finite jump in 
(e) at Tn = 124 K. An investigation of the probability distribution (not shown) of the internal energy 
at Tn exhibits the two maxima structure that is typical for first order phase transitions. 

For D = J, we find a sudden collapse of the magnetic moment at Tg = 127 K which is connected 
with an abrupt decrease of the volume of about 3 %. The magnetic phase transition takes place a few 
Kelvin below at Tn = 120 K. This situation is reflected in a double jump of the internal energy. Upon 
cooling down, a large thermal hysteresis is encountered: The system has to be cooled down to — 53K 
before it jumps back into the AF HM state. Responsible for this hysteresis is that lattice strain imposes 
a large free energy barrier between HM and LM states which is unlikely to be surmounted by our 
constant pressure Monte Carlo algorithm. However, the large energy barrier is suppressed in Fig. ^ 
since at constant volume the system is forced into a mixture of two coexisting HM and LM phases. 
Similar first order transitions with a large thermal hysteresis have also been observed in simulations of 
ferromagnetic systems; for a detailed discussion we refer to |6|. Above Tg the mean magnetic moment 
is considerably lower than the high temperature limit {S^') = 2/3 of a spin-1 Ising model. So with 
increasing temperature more atoms will be excited to the HM state, leading to an enhanced thermal 
expansion. This process is further encouraged with increasing temperature, because the displacement of 
the atoms due to lattice vibrations and the difference between HM and LM neighbor distances become 
comparable. This levels out the energetic difference between Si = and Si = ±1 states introduced by 
Hv For the same reason, discontinuous magnetovolume effects do not occur for smaller values oi D: At 
higher temperatures, the HM and LM minima of the free energy are smeared out by lattice vibrations. 

Above D = 2 J, the ground state of the system is nonmagnetic. So HM-LM transitions as well 
as magnetic phase transitions do not appear. However, we observe an enhanced thermal expansion 
coefRcient, due to excitation of HM atoms. 

DISCUSSION 

For D = J, our model shows magnetovolume effects similar to the properties of YMn2 mentioned in the 
introduction: An abrupt collapse of the mean magnetic moment connected with a decrease of the volume 
of several percent and a largely enhanced thermal expansion coefficient. Lowering D produces continuous 
magnetovolume effects. This resembles the situation in Y(Mni_a;Ala;)2, where small concentrations 
of Al destroy the first order transition and anti- Invar behavior |l2|. With increasing Al content the 
magnetovolume effects gradually vanish. This is due to the enlargement of the lattice constant by the 
Al atoms which makes it energetically more difficult for the Mn atoms to switch to the LM state. On 
the other hand, substituting a few percent of the Y atoms with the smaller Sc has the opposite effect: 
The Mn atoms are nonmagnetic at T = 0, the Mn moment increases with temperature, leading to an 
enhanced thermal expansion coefficient . This is the situation we observed for Z) > 2 J. 
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Figure 1: Energy per lattice site and average magnetic moment at T = 4 K as a function of the lattice 
constant. 
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Figure 2: Energy per lattice site for various values of D as a function of the temperature. 
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Figure 3: Temperature dependence of the average magnetic moment for various values of D. 
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Figure 4: Temperature dependence of the volume and the hnear thermal expansion coefficient for various 
values of D. 
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